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ABSTRACT 

Context. The winds of massive stars create large (> 10 pc) bubbles around their progenitors. As these bubbles expand 
they encounter the interstellar coherent magnetic field which, depending on its strength, can influence the shape of the 
bubble. 

Aims. We wish to investigate if, and how much, the interstellar magnetic field can contribute to the shape of an 
expanding circumstellar bubble around a massive star. 

Methods. We use the MPI-AMRVAC code to make magneto-hydrodynamical simulations of bubbles, using a single star 
model, combined with several different field strengths: B = 5, 10, and 20 pG for the interstellar magnetic field. This 
covers the typical field strengths of the interstellar magnetic fields found in the galactic disk and bulge. Furthermore, 
we present two simulations that include both a 5 pG interstellar magnetic field and a warm (10 000 K) interstellar 
medium (ISM) and two different ISM densities to demonstrate how the magnetic field can combine with other external 
factors to influnece the morphology of the circumstellar bubbles. 

Results. Our results show that low magnetic fields, as found in the galactic disk, inhibit the growth of the circumstellar 
bubbles in the direction perpendicular to the field. As a result, the bubbles become ovoid, rather than spherical. Strong 
interstellar fields, such as observed for the galactic bulge, can completely stop the expansion of the bubble in the 
direction perpendicular to the field, leading to the formation of a tube-like bubble. When combined with an ISM that 
is both warm and high density the bubble is greatly reduced in size, causing a dramatic change in the evolution of 
temporary features inside the bubble such as Wolf-Rayet ring nebulae. 

Conclusions. The magnetic field of the interstellar medium can affect the shape of circumstellar bubbles. This effect 
may have consequences for the shape and evolution of circumstellar nebulae and supernova remnants, which are formed 
within the main wind-blown bubble. 

Key words. Magnetohydrodynamics (MHD) - Stars: circumstellar matter - Stars: massive - ISM: bubbles - ISM: 
magnetic fields - ISM: structure 


1. Introduction 

Massive stars lose a significant fraction of their mass in 
the form of stellar wind, which expands into the surround¬ 
ing medium. As it expands, the stellar wind collides with 
the gas in the interstellar medium (ISM), creating a low 
density bubble, expanding over time. Analytically, the gen- 
eral shape of s uch a wind blown bubble was p redicted by 


ISM and time-dependent wind parameters (E.g. 

jRozyczka 

1985; Rozyczka & Tenorio-Tagle 

1985b ja|). Eurtl 

ler devel- 


lution theory made it possible to follow the evolution of the 
stellar wind parameters as the progenitor star evolves, and 
to investigate the effects on the morphology of circumstellar 
bubbles. 


Avedisova (e.g 1972) and Weaver et al 


VVCIO pj. uy ;- - 1 - 

( 129771 ) However Initial models by [Garcia-Segura et al.| ( |1996b|a[ ) fo- 


these models were limited, because they were strictly one- cussed primarily on the transitional phases of the stellar 


dimensional and could therefore not take into account the 
effects of either a-spherical stellar winds, or irregularities 
in the ISM. Furthermore, they assumed that the stellar 
wind properties would remain constant over long periods of 
time. These analytical models were complemented by nu¬ 
merical simulations, which can include both a non-uniform 
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evolution, e.g. the transition from main sequence to either a 
Luminous Blue Variable or a Red Supergiant and the tran¬ 
sition from the giant stage to the Wolf-Rayet stage. These 
models proved that the radical changes in stellar wind pa¬ 
rameters accompanying such transitions lead to the forma¬ 
tion of shells that can be observed as circumstellar nebulae. 
Subsequent models by 


Marie et al. (2005 2007 


Freyer et aL| ( |2003[ |2006[ ) and 


van 


2008) followed the entire evolution 
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of the star and introduced the e ffect of pho t o-ioni z ation . 


Tenorio-Tagle et ah (1990) and Dwarkadas (2005, 2007) 
added the supernova explosion and followed tne evolution 
of the su pernova remnant insid e the circumstellar bubble. 
Recently, Toala & Arthur (2011) included radiative transfer 
to investigate the ionization structure of the circumstellar 
medium and added thermal conduction to their models. 


Geen et al. (2015) made similar models for a somewhat 
lower mass star (15 Mq) in 3D. 

Another factor potentially affecting the evolution of 
circumstellar bubbles is the interstellar magnetic field. 
Observations show that the interstellar medium contains 
large scale coherent magnetic fields. These fields ar e usn- 
ally rather weak {1 < B < 60 pG according to e.g. |Rand 


fc Knlkarni||1989| Ohno fc Shibata|1993| [Frick et al.|2001 


Opher et al.||2Q09 ; Shabala et al.||2Q10[ [Crocker et al.||2011 


Fletcher et al. |20li||Heerikhuisen fc Pogorelov|2011||Vallee| 


201 1| [Green et al.||2bl2| |, but can potentially influence the 

shape of the interstellar bubbles, blown by stellar winds. 
(Higher values have been found for the interior of molecular 
clouds, though these typically occur in those regions where 
the cloud is at its densest, r e nderin g them relatively small 


in scale. Even so, Crutcher (1999) shows a field strength 
of 480 pG with a radius of 22 pc for SgrB2.) Despite their 
low strength, interstellar magnetic fields can influence the 
shape of circumstellar bubbles. This influence can occur on 
a relatively small scale, since the presence of a magnetic 
field changes the properties of hydro dynami cal instabili- 


ties (lJun et al. 1995 [Stone fc Gardiner [200 7). This effect 
was aescrib ed analytically ( jPgani et^ JM [Breitschwerdt| 
et al. 2000) and shown in numerical models of a stellar wind 


bow shock by van Marie et al. (2014b). On a larger scale, the 
interstellar niagnetic held can influence the general shape 
of a circumstellar bubble by reducing the expansion rate 
in the direction perpend icular to the mag netic field. This 
was already predicted by Heiligman (1980) and wo rked out 
in detail for the evolution of planetary nebulae byjSoker fc 



De ^ani| ( |T9^ . More recently, [Falceta-Goncalves fc Monteiro| 


showed numerically that a strong interstellar mag- 
netic field 500 pG) could lead to the creation of a bipolar 
planetary neb ula. A similar effect was noted by [van Marlej 
et al. ( 2014a| ), who used numerical simulations to demon¬ 
strate that interstellar magnetic fields can explain the for¬ 
mation of “eye-like” struct ures around cert ain asymptotic 
giant branch (AGB) stars ( Gox et ar]|2012 ). 

Hot, massive stars have strong winds, which might allow 
them to overcome the forces exerted by the magnetic field. 
On a large scale, the influence of the interstellar magnetic 
held on the shape of super bubbles, which conta in clusters 
of massive st ars, was investigated nume rically by [Tomisakaj 


(1990, [1992) and Ferriere et al. (1991) who found that a 


parallel magnetic field constrains the expansion o f the bub¬ 
ble in the direction perpendicular to the held. [Tomisak^ 
(1998) investigated the 3D dynamics of a superbubble in 


stratified and magnetized medium, exploring the conditions 
under which the bubble can blow out in the halo. However, 
these models considered supernova explosions to be the only 
source of energy in the bubble and approximated the effect 
of these sequential supernovae with an input of energy that 
remained constant in time. No attempt was made to in¬ 
clude the stellar wind as a separate force to drive the initial 
expansion of the bubble, or to investigate how the devia¬ 
tion from spherical symmetry of the outer shell influences 
the morphology of the circumstellar shells formed inside 


the bubble when the properties of the stellar wind change 
or during the expansion of a supernova. Although the to¬ 
tal energy contribution of the supernova to the circum¬ 
stellar bubble is approximately three tim es larger than the 
wind energy contribution (See Table 2.2), the wind shapes 
the bubble in which the supernova expands and therefore 
should not be neglected in the simulations. We use the MPI- 


AMR VAC magneto-hydrodynamics code (Keppens et al. 


2012) to run a series of simulations of circumstellar bub¬ 


bles, using a generic 40 Mq star model like the one used ir 


van Marie et al. 


( 2012 ^ 


By varying the strength of the in- 
terstellar magnetic field, while keeping all other quantities 
constant, we investigate both its qualitative and quantita¬ 
tive effects on the wind expansion to see if, and how, the 
magnetic held influences the shape and evolution of the cir¬ 
cumstellar bubble. Additionally, we present two simulations 
showing similar expansions in a warm, rather than a cold 
ISM to demonstrate how the interstellar held can combine 
with other external factors to influence the size and shape 
of circumstellar bubbles. Instead of using a constant inflow 
of en ergy at the centre, as was done by [Tomisakaj ([1990 
1992) and Ferriere et al.[ ( [1991 ), we approximate the stel¬ 
lar evolution with a three-stage model (main sequence, red 
supergiant, and Wolf-Rayet) followed by a supernova explo¬ 
sion. This allows us to follow the evolution of the temporary 
shells formed inside the bubble as a result of changes in the 
wind parameters and the eventual supernova explosion. 


1.1. Layout 

In Sect. of this paper we describe the physical consid¬ 
erations that are involved in the interaction between the 
expanding stellar wind and the interstellar magnetic held 
and use an analytical approximation to predict the field- 
strength necessary to overcome the ram-pressure of the 
wind. We then describe our model parameters and numer¬ 
ical method in Sect. [U The results of our simulations for a 
cold ISM are shown in Sect. [4j The case of a warm inter¬ 
stellar medium is shown in Sect. El All cases are discussed 
in Sect.j^ Finally, in Sect.[^ we present our conclusions. 

For full animations of our results in electronic form, see 
Appendix 


2. The interaction between a wind-blown bubble 
and the interstellar magnetic field 

2.1. Hydrodynamical interaction 


The physics that govern the expansion of a wind-blown bub¬ 
ble are well understood. As the stellar wind collides with the 
ISM its kinetic energy is converted to thermal energy and 
it forms a bubble of hot, shocked gas. This bubble is con¬ 
strained on the inside by the wind termination shock, where 
the ram-pressure of the wind (pram) balances against the 
thermal pressure in the hot bubble (px), and on the outside 
by a shell of compressed interstellar gas that i s pushed out¬ 


ward s by the thermal pressure of the bubble (Weaver et al. 
1977 ). (N.B. In this section, all pressures related to the 
ISM are denoted by an upper case P, whereas all pressures 
related to the wind are denoted by lower case p.) As the 
shell expands, its motion is determined by the balance be¬ 
tween three different pressure sources. On the inside of the 
shell is the thermal pressure of the hot, shocked stellar wind 
bubble. On the outside the shell is confined by the thermal 
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Fig. 1. Schematic of a wind-blown bubble inside a galac¬ 
tic magnetic field. As the bubble expands it compresses 
the magnetic field lines (black lines), which increases the 
local magnetic pressure. The expansion of the outer shell 
(solid red line) stops when the thermal pressure (px), which 
drives the expansion is equal to the external magnetic field 
pressure Pb- This determines the location of the wind ter¬ 
mination shock (dashed red line), which is located at the 
point where the ram pressure of the wind (pram) is equal 
to pt- 


pressure of the ISM and the ram-pressure of its own motion 
into the ISM. Of the two confining pressures, the thermal 
pressure of the ISM is generally considered negligible com¬ 
pared to the ram-pressure. This is based on the assumption 
of a cold ISM (temperature ^ 100 K). The thermal pres¬ 
sure Pt is determined by the ISM density (pism) and the 
temperature (Tjsm) as 


PT = ^^kTiSM, ( 1 ) 

mhPe 


with mh the hydrogen mass, pe the mean atomic weight 
per free electron and k the Boltzmann constant. The ram- 
pressure experienced by the shell in the co-moving frame is 
equal to 

-Pram = PlSM'^shelb (^) 


with 'T’sheii the expansion speed of the shell. Assnr ning an ex¬ 
pansion velo city of 'i^sheii of about 10 kms“^ (e.g. van Marie 


et al. 2012), the ram-pressure is two orders of magnitude 
above the thermal pressure. 

In a warm ISM (8 000-10 000 K) the situation is differ¬ 
ent, since the thermal pressure exerted by the ISM is of 
the same order as the ram pressure. Furthermore, assuming 
that the swept-up shell will not cool below the temperature 
of unshocked ISM, the shell will not be compressed. Under 
those circumstances, it is quite possible that the outward 
motion of the expanding bubble will become subsonic. If 
this occurs the forward shock will disappear and the shell 
will become an outward-moving eddy, with only a slight 
density increase over the ambient ISM. The contact con¬ 
tact discontinuity will remain, but rather than separating 


shocked wind from shocked ISM, its location will be deter¬ 
mined by the pressure balance between the thermal pres¬ 
sure of the shocked wind and the thermal pressure of the 
unshocked ISM. A further complication arises if the warm 
ISM is in fact heated by the radiation from the star, which 
can form an HII region outside its win d-blown bubble as 
demonstrated in the numerical models of|F reyer et al. ( 20031 

Tbala & Arthur 


2QQ6|); van Marie et al. (2005 2007) and 


]2011). Such an w region expands outward into the cold 


ISM, even as it is being swept-up by the stellar wind, giv¬ 
ing rise to a layered bubble consisting partially of shocked 
wind and partially of photo-ionized ISM. 


2.2. Effect of an interstellar magnetic field 

The mag netic fields i n the ambient medium can vary con¬ 
siderably (Beck 2QQ9|) . In the galactic disk, field strengths of 
5-10 pG seem to be the norm, whereas in the galactic bulge, 
the strength of the field can be more than twice as high 
( Rand fc Knlkarni 1989 Ohno fc Shibata|[l993[ [Shabala 
et al.|2010 ). Although these helds are not very strong in an 
absolute sense, they extend over large scales (10-100pc). If 
the ISM contains a magnetic field, the expanding hot bub¬ 
ble has to overcome a constraining force comprising three 
separate components: 

1. The thermal pressure of the ISM 

2. The ram pressure created by its own expansion into the 
ISM 

3. The magnetic field pressure of the ISM (for the expan¬ 
sion component perpendicular to the magnetic field). 

Furthermore, the nature of the outward expansion is no 
longer determined exclusively by the sound speed in the 
ISM. Instead, the expansion will manifest itself as a fast 
magnetosonic wave that travels into the ISM. As long as the 
wave is both supersonic and super-Alfvenic, the swept-up 
ISM will be compressed, resembling the shells predicted by 
the purely hydro dynamical model. In a cold ISM the Alfven 
speed, even for a 5pG magnetic field is higher than the 
speed of sound. Hence, the expansion wave will become sub- 
Alfvenic, even while still supersonic. In terms of the shock 
conditions, this means that the shock makes a transition 
from a J-type shock to a G-type shock, because the fact 
that the Alfven speed exceeds the speed of sound allows 
the Alfven waves to carry information ahead of the shock. 
Once this happens, the compression will be reduced. 

Assuming that the expansion is (at least initially) super- 
Alfvenic, the pressure in the compressed shell will be deter¬ 
mined by a combination of thermal pressure due to shock 
heating and, depending on the angle of the expansion with 
the direction of the magnetic field, the magnetic pressure 
of the compressed magnetic field. If the magnetic fiux that 
originally passed through a surface area 4x7?^, with R the 
outer edge of the expanding shell of shocked, ISM is com¬ 
pressed into a shell with thickness the magnetic field 
strength in the shell (Hgheii) will become 


-Rshell = 




ISM 


47r(R2 _ {R-DY) 


= -Rism 


R^ 


2RD - ’ 


(3) 


in the direction perpendicular to the field lines. Because 
magnetic pressure scales with the field strength squared, a 
high compression rate {D/R 1/10), would indicate an 
increase in magnetic pressure of ^ 25. This increase in 
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magnetic pressure has to be balanced at the inner bound¬ 
ary of the shell by a higher thermal pressure in the shocked 
wind, which in turn requires the wind termination shock to 
be closer to the star. This also means that the expansion 
velocity of the hot bubble has to be reduced (the shocked 
wind decompresses to drive the expansion). At the outer 
edge of the shell, the magnetic pressure can only be bal¬ 
anced by an increase in the ram pressure undergone by the 
expanding shell, which would require an increase in expan¬ 
sion velocity at the forward shock. As a result, the swept-up 
shell starts to decompress. The amount of decompression 
depends on the local strength of the magnetic pressure, 
which in turn depends on the angle of the shell’s motion 
with the magnetic field. Consequently, we can expect the 
contact discontinuity to follow an ovoid shape, with the ma¬ 
jor axis along the direction of the magnetic field. The shell 
will have a varying thickness, being thickest along the mi¬ 
nor axis of the ovoid. Because of the decompression of the 
shell, the difference in magnetic field pressure between the 
shell and the ambient medium is reduced. Eventually, the 
shell, which by this time will have become sub-Alfvenic, will 
decompress so far, that its magnetic field pressure is effec¬ 
tively equal to the ambient magnetic field pressure. At that 
moment we can no longer speak of a shell. Instead, the hot, 
shocked-wind bubble will be in direct pressure equilibrium 
with the ambient medium and the expansion perpendicu¬ 
lar to the magnetic field stops. This will eventually occur 
for any wind-driven bubble expanding in a magnetic ISM. 
Obviously, the question is whether the star will live long 
enough for its bubble to reach this point. 

In order to estimate whether the circumstellar bubble 
will reach the point where it stops expanding perpendic¬ 
ular to the magnetic field and at what distance this will 
occur, we compare the two opposing forces that must reach 
an equilibrium: the ram pressure of the stellar wind, and 
the magnetic field pressure of the ISM. Albeit that the ram 
pressure of the stellar wind near the surface of a massive 
star is much stronger than the interstellar magnetic field, it 
is proportional to the density (p), which in a free-streaming 
wind decreases with IjR^. As the stellar wind expands, it 
approaches the point where its ram pressure (pram) equals 
the interstellar magnetic field pressure ( Pr)- We can deter¬ 
mine the distance at which this takes place: 


Table 1. Wind and supernova parameters, with tend the 
time from the start of the simulation to the end of each 
evolutionary phase, M the wind mass loss rate, the wind 
velocity, and M and E the total mass and energy injected 
into the medium. 


Phase 

tend 

M 

Vw 

M 

E 


[Myr] 

[ Mo yr-i] 

[kms“^] 

[Mo] 

10"" [erg] 

40 M© 

MS 

4.3 

1.0x10“® 

2 000 

4.3 

1710 

RSG 

4.5 

5.0x10“® 

15 

10.0 

0.224 

WR 

4.8 

1.0x10“® 

2 000 

2.7 

1190 

SN 




10.0 

10 000 


wind (M = 10“^ Mq yr“^ and Vo© = 2000 km s“^) and a 
magnetic field strength of bpG, we find that Rb — 7.3 pc. 
This is rather large. However, massive stars may produce a 
bubble large enough to get to this point, depending on the 
density of the ISM. (Eldridge et al. 2006 found typical re¬ 
verse shock radii of ^ 5 pc for massive stars in a 10 # cm“^ 
ISM at the end of their evolution.) Therefore, although the 
magnetic field pressure will influence the shape of the bub¬ 
ble, it will in most cases not stop the expansion. On the 
other hand, at P = 20 pG, Rb — 2.6 pc, placing it within 
the termination shock of a typical 0-star or Wolf-Rayet star 
wind. 

Beside the magnetic pressure, we also have to consider 
the magnetic tension force. As long as the field lines are 
parallel, the magnetic tension force is zero. However, when 
the circumstellar bubble expands, it will exert a force on 
the field lines, forcing them into a curved shape. The mag¬ 
netic tension force, which is directed inward with respect 
to the curvature of the field lines, will try to counteract the 
expansion of the bubble. Therefore, the expansion of the 
bubble will always be reduced in the direction perpendicu¬ 
lar to the field, even if the magnetic pressure is too small 
to overcome the ram pressure of the wind. 


3. Method 

3.1. Computational method 


Pram 

= Pb, 

(4) 



(5) 


2/io ’ 

M Voo 


(6) 



47rR| 


Rb 

= ^\/2Mvoo, 

(7) 

with Voo 

the terminal wind velocity, B = 

B the magnetic 


field strength po = Air the magnetic permeability of vac¬ 
uum in cgs. units, M the mass loss rate of the star and 
Rb the radius at which the ram pressure of the stellar wind 
equals the magnetic field pressure in the ISM. Once the 
wind termination shock reaches this distance from the star, 
the ram pressure of the stellar wind (and therefore the ther¬ 
mal pressure of the shocked wind) becomes less than the 
interstellar magnetic field pressure. As a result, the bubble 
will stop to expand in the direction perpendicular to the 
magnetic field (see Eig. [^. Eor an 0-star main sequence 


We use the MPI-AMRVAC hydrodynamics code (Keppens 


et al. 2012), which solves the conservation equations for 

0 , 


mass. 


1 + V^(,v) 


(8) 


momentum: 


1 


Vptot-(V X B) X B = 0, (9) 

Mo 


and energy: 

l + V-(ev) + V.(p„v) = -(;^)V), (10) 


with V the velocity vector, ptot fhe sum of the thermal 
pressure and the magnetic pressure, and the hydrogen 
mass. The energy density e is defined as the sum of thermal, 
kinetic and magnetic energy density. 


e 


P 

7-1 



2/io ’ 


(11) 
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with 7 = 5/3 the adiabatic index and v = | v|. We assume 
ideal MHD, so there are no local source terms for the mag¬ 
netic field. Consequently, 


^ - Vx (vx B) = 0. 


( 12 ) 


We solve these equations using the total variation dimin¬ 
ishing, Lax- Friedrich TVDL F method combined with a koren 
flux limiter ( Kuzmin|2QQ6 ), except for the two highest levels 
of refinemen t, where w e apply the more diffusive minmod 
fiux limiter ( Roe||1986 ) to prevent the growth of numerical 
instabilities at the polar axis. 

The divergence of the magnetic field is kep t at zero with 
the method described by Powell et al. (1999). Equation [l^ 
includes the effect of radiative cooling, which depends on 
local density, as well as a temperature dependent cooling 
curve for sol ar-metallicity gas, A(T), obtained from Schure 
et al. (2009). This curve, generated with the SPEX code, 
covers a temperature range of 10^-10^ K. For higher temper¬ 
atures, we assume that the cooling scales with because 
the Bremsstrahlung dominates the cooling. For tempera¬ 
tures below T = 10^ K, we copy the approach of [Schure 


et al. 


_ _ ([20091), by using the cooling function from [Dalgarno 

&: McCray ' ^1972 ), adjusted for a pre-set ionization fraction, 
which we choose to be 10“^. We set a minimum tempera¬ 
ture of 100 K throughout the simulation to prevent numeri¬ 
cal problems that can arise because of extreme compression 
of radiatively cooling gas. 

The simulations are set up in a manner similar to the 


one used in van Marie et al. (2012), with a 2D cylindrical 
grid, which is symmetric around the z-axis. We start with 
a basic resolution of 0.35 pc per grid cell and allow 3 ad¬ 
ditional levels of refinement, each doubling the resolution. 
This gives us an effective resolution of 0.044 pc. This grid 
is filled with ambient medium according to the parameters 
specified in Sect. 3.2 The stellar wind is introduced in the 
grid by filling a small half-sphere (0.25 pc in radius) on the 
z-axis with wind material, with the wind parameters chang¬ 
ing over time according to the evolutionary phase of the star 
(Table 1). Around this sphere we enforce an additional level 
of grid refinement (0.022 pc per cell) to reduce numerical er¬ 
rors, a possible occurrence, when a spherical expansion is 
simulated on a rectangular grid. The stellar wind parame¬ 
ters are the s ame as those used for t he 40 Mq star in our 
earlier paper (van Marie et al. 2012). The mass loss r ates 
and evolution tim e are based on the model used in IvanI 
Marie et al.| ( |200'5 ). Wind velocities reflect those typically 


found in observations (Earners & Cassinelli 1999, and ref¬ 


erences therein). The mass-loss rate in the final phase (the 
Wolf-Rayet phase) has been reduced to refl ect the lower 
values obtained by [Vink fc de Koter (2005) when taking 
into account the effects of clumping. 

The supernova is introduced in a similar manner: during 
the first time step of the supernova simulation, the same 
sphere is filled with hot gas with a total mass of 10 Mq 
and a total thermal energy of lO^^erg. The supernova mass 
reflects the final mass of the star as well as the values cal¬ 


culated by, e.g., Eldridge fc Tout (2004). For simplicity, we 
use a supernova energy o f 10^^ei^, which is typi cal of core¬ 
collapse supernovae (e.g. van Marie et al. 2010[ and refer¬ 
ences therein). During the supernova phase we allow the 
maximum level of resolution (0.022 pc per cell) throughout 
the entire domain to make sure we resolve the high velocity 
shock front. 


Throughout the simulations, the interstellar magnetic 
held Bism is set parallel with the z-axis and fixed at the 


outer boundari es, using the same method as in [van Marie 
al. ( 2014b|a ). We do not tak e into consideration a stellar 
mag netic held. Simulations by ^d-Doula fc Owocki (2002) 
and |Ud-Doula et al.| ([20081 |2009p show that the magnetic 
fields of massive stars, important when close to the stellar 
surface, are typically dominated by the ram pressure of the 
wind at large distances. 


3.2. ISM parameters 

We run four simulations varying only the magnetic held 
strength in the ISM: 

— Simulation A, the “basic” model. A 40 Mq star in an 

ambient medium of 10“^^‘^gcm“^ 20#cm“^ with¬ 

out a magnetic held. The ambient medium temperature 
is set to 100 K. This model serves as a reference against 
which we measure the other models. 

— Simulation B, similar to simulation A, but with an in¬ 
terstellar magnetic held of bpG, which is representative 
of the weak magnetic fields in the galactic disk. 

— Simulation C, similar to simulation B, but with an inter¬ 
stellar magnetic held strength of 10 pG, representative 
of the strong magnetic fields in the galactic disk and the 
weak fields in the galactic bulge. 

— Simulation D, similar to simulations B and G, but with 
an interstellar magnetic held strength of 20 pG, repre¬ 
sentative of the strong magnetic fields in the galactic 
bulge. 

These allow us to explore both the qualitative and quan¬ 
titative effects of the influence of the interstellar magnetic 
held on the morphology of the circumstellar bubble, and, 
additonally, we run two simulations for an ISM with a tem¬ 
perature of 10 000 K to show how the effect of the interstel¬ 
lar magnetic held combines with other forces acting on the 
bubble: 

— Simulation E, similar to simulation B but with an inter¬ 
stellar density of 10“^^‘^gcm“^ and an ISM tempera¬ 
ture of 10 000 K as can be expected in warm, low density 

ISM. 

— Simulation F, similar to simulation B, but an ISM tem¬ 
perature of 10 000 K, reminiscent of the conditions in an 
Hii region. 


4. Bubbles in a cold ISM 

4.1. Expansion in the absence of a magnetic field 


Simulation A shows our basic model, which was also used in 

If 


van Marie et al 
bubble follows t 


(2012). The evolution of the circumstellar 


e expected path. The main sequence wind 
creates a large, spherical bubble of hot, shocked wind ma¬ 
terial, which expands into the ambient medium, sweeping 
up a sh ell, following the pa ttern predicted by the analytical 
model ([Weaver et al. |l977l ). This is shown in the left panel 
" The shell is subj ect to instabiliti es, primarily of 

shell type (Vishniac 1983), in the earlier 


of Fig 

the linear thin 
stages. These instabilities result from the different forces 
acting on the shell: isotropic thermal pressure on the inside 
and ram-pressure from the outside. Rayleigh-Taylor insta¬ 
bilities, resulting from the density difference at the contact 
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Fig. 2. Logarithm of the density [cgs.] for the circumstellar medium in the absence of an interstellar magnetic field. The 
left panel shows the circumstellar bubble on a 40x40 pc frame after 1 (left side) and 2 (right side) Myr have past since 
the start of the main sequence. The bubble is spherically symmetric, except for local thin-shell and Rayleigh-Taylor 
instabilities. In the centre (frame size 60x60pc), the star leaves the main sequence after 4.3 Myr (left side) and becomes 
a red supergiant, which lasts until 4.5 Myr after the start of the simulation (right side), forming a new shell at the wind 
termination shock. The right panel (frame size 70x70 pc) shows the transition to the Wolf-Rayet stage, which forms a 
new shell that sweeps up the red supergiant wind (left panel, lOOOOyrs into the Wolf-Rayet stage) and collides with the 
red supergiant shell. The fragments of this collision move out into the old main sequence bubble (right side, 100 000 years 
into the Wolf-Rayet fase). 



R-axis (xlO^lS cm) R-axis (xlO^lS cm) R-axis (xlO^lS cm) 

Fig. 3. Similar to Fig. showing the expansion of the supernova remnant. All frames have the same size of 80 x 80 pc. 
At the end of the Wolf-Rayet phase (left side of left panel) the density in the hot bubble has become smooth once again. 
The supernova starts as a high pressure blast wave at the centre of the bubble and quickly expands. After 400 years it 
has swept up nearly half of the free-streaming Wolf-Rayet wind (right side of left panel). 4000 years after the explosion 
(left side of centre panel), the supernova remnant has swept up most of the hot bubble. After 7200 years (right side, 
centre panel) it collides with the outer shell. Afterwards, the supernova remnant flows back toward the centre of the 
bubble (12 000 years after the explosion, left side of right panel). Eventually, its motion will slow down to subsonic level 
and the pressure in the bubble will start to equalize (right side of right panel, 40 000 years after the explosion). 


discontinuity between the shocked wind and the swept-up 
shell develop more slowly, but begin to dominate in the later 
stages of the main sequence (right side of the left panel and 
left side of the centre panel of Fig. |^. 

As the star makes the transition to the RSG phase, the 
wind parameters change, leading to a change in the inner 
part of the CSM. The high density RSG wind forms a shell 
at the wind termination shock (right side of the centre panel 
of Fig. [^. Eventually, the star becomes a Wolf-Rayet star, 
with a mst, powerful wind that sweeps up the RSG wind 
material in a shell (left side of the right panel of Fig. [^. 


This Wolf-Rayet wind driven shell collides with the RSG 
wind shell and both shells are fragmented. The individual 
fragments move out into the shocked main sequence wind 
bubble, where they dissipate or merge with the outer shell 
(right side of the right panel of Fig. 


Segura et al. 

(e.g. 1996a); van Marie et al. (e.g. 

20051); 

Freyer et al.| ( 

e.g. 2006); Dwarkadas (e.g. 2005) and Toala 
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t [Myr]' 


Fig. 4. Values for the wind termination shock obtained 
from simulation A (black) at the end of each evolution¬ 
ary phase, compared to Rb for the three magnetic field 
strengths. At SpG (red), the magnetic field can never stop 
the expansion of the circumstellar bubble perpendicular to 
the field and the 10 pG (green) field will do so, but only at 
the very end of the main sequence. The 20 pG (blue) field 
will stop expansion throughout most of the stellar evolu¬ 
tion. 


Using the results from this simulation we can make a 
rough prediction as to what extent the magnetic fields in 
the subsequent simulations will inhibit the expansion of the 
circumstellar bubble. Figure shows the wind termination 
shock at the end of each evolutionary phase, compared to 
the Rb values for the stellar wind parameters and various 
magnetic field strengths. As long as the termination shock 
radius is smaller than Rb the wind driven bubble will be 
able to expand perpendicularly to the magnetic field, albeit 
at a lower speed owing to the magnetic tension force and 
the increased pressure from the ISM. This is what we can 
expect for simulations B and G. For simulation D, which has 
an Rb that is lower than the termination shock for most 
of the simulation, expansion perpendicular to the field will 
stop completely. 

Once the star reaches the end of its evolution, it ex¬ 
plodes as a core-collapse supernova (left panel of Fig. [^. 
This causes a sudden injection of mass and energy at the 
centre of the circumstellar bubble, which by now has a ra¬ 
dius of approximately 30 pc. The supernova remnant ex¬ 
pands quickly into the low density medium until it reaches 
the inner edge of the swept-up ISM shell (centre panel of 
Fig.l^. There it is stopped by the high mass (^ 40 000 Mq) 
sheT Failing to break out, the supernova remnant recoils 
and starts an inward motion (left side of right panel of 
Fig. [^, which continues until it reaches the centre of the 
bublne. From this point the movement slows down until, 
finally, the density and pressure inside the bubble start to 
even out (right side of right panel of Fig. [^. 


tive of the magnetic field i n the spiral arms of the galac- 


tic disk on 10-10 0 pc scales ([Rand fc Kulkarni|| 19891 
& Shibata 1993 Shabala et al. ' 2 QTq) ~ Green et ah' 


Ohno 

2 ^ 


and the l ower values obtained for magnet ic helds in nearby 
galaxies ( [Beck|20Q9 Fletcher et al.|2QlT ). During the main 
sequence phase, the bubble expands outward, creating a 
non-magnetic hole in the field by pushing the field lines, 
embedded in the ISM, outwards (left and centre panels of 
Figj^. The magnetic tension force, which is directed inward 
with respect to the curvature, counteracts the displacement 
of the field lines, reducing the expansion perpendicular to 
the field lines. At the same time the field lines are pressed 
together in the shocked shell, increasing the magnetic pres¬ 
sure, but its strength is insufficient to stop the expansion. 
This results in an ellipsoid circumstellar bubble. 

Owing to the asymmetrical shape of the bubble, the 
remnants of the WR-RSG shell collision will first encounter 
the outer shell at the point where it has expanded least (the 
direction perpendicular to the external field). This causes 
the shell fragments to be pushed aside along the z-direction, 
causing a mirror symmetry, but most of the turbulence dis¬ 
appears before the star reaches the end of its evolution 
(right panel of Fig. and left side of left panel of Fig. 

At the end of the stellar evolution (left side of left panM 
of Fig. the shocked wind bubble is elliptical, with an 
asymmetry of approximately 1.3:1, because the magnetic 
field pressure limits the expansion in the direction perpen¬ 
dicular to the field (r-direction). Along the axis parallel 
with the magnetic field, the bubble has a radius of approx¬ 
imately, 30 pc, comparable with the radius of the bubble in 
simulation A, whereas the radius along the shorter axis is 
approximately 23 pc. 

The magnetic field lines that are swept up in the outer 
shell resist being pressed together, consequently the com¬ 
pression in the shell is lower than in the model without 
a magnetic field. This results in a shell with a variable 
thickness: widest where the shell moves perpendicular to 
the field and thinnest where it moves parallel to the field, 
resulting in the outer edge of the shell being nearly spher¬ 
ical, despite the elliptical shape of the contact discontinu¬ 
ity. The shell lacks instabilities, because the magnetic pres¬ 
sure makes it too thick for thin-shell instabilities (Fig[^. 
Rayleigh-Taylor instabilities could still develop, but they 
would have to move perpendicular to the magnetic field 
lines. Such instabilities can be seen faintly along the con¬ 
tact discontinuity (best observed in centre panel of Fig. [^, 
but the instabilities do not grow over time. 

Eventually, the star explodes as a supernova. The evo¬ 
lution of the supernova remnant (Fig. closely resembles 
that of the non-magnetic models, since there is no direct 
contact between the magnetic field and the supernova rem¬ 
nant, which is contained by the shocked wind bubble. Like 
the expanding WR shell, the supernova first encounters the 
outer shell at the “waist” of the bubble (right side of cen¬ 
tre panel in Fig. 1^, causing a mirror symmetric structure 
inside the hot bubble. Over time this structure tends to 
disappear until, eventually, the density inside the bubble 
becomes nearly constant (right panel of Fig. [^. 


4.2. Low interstellar magnetic field: 5luG 

Simulation B shows the evolution of the circumstellar bub¬ 
ble in a weak (bpG) external magnetic field, representa¬ 


4.3. Intermediate interstellar magnetic field: lOjuG 

During the main sequence evolution. Simulation G, which 
uses a magnetic field of 10 pG (high for t he spiral arms, but 


low for the galactic bulge according to Rand fc Kulkarni 
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Fig. 5. Similar to Fig. but for Simulation B with an external magnetic field of 5[JiG, showing the progression of the 
bubble evolution at 1, 2, 4.3, 4.5, 4.51 and 4.6 Myr. Unlike the non-magnetic simulation, the swept-up shell of shocked 
ISM is relatively thick, which prevents the formation of thin-shell instabilities. The overall-shape of the bubble is ellipsoid 
owing to the magnetic pressure which slows down expansion perpendicular to the field, but the evolution of temporary 
structures inside the bubble is similar to the non-magnetic case. A faint sign of instability can be seen along the contact 
discontinuity, but the instabilities do not grow over time. 



Fig. 6. Similar to Fig.^ but for Simulation B with an external magnetic field of bpG. Initially, the supernova expansion 
is nearly identical to the non-magnetic case, but since the supernova expansion is spherical and the outer boundary of 
the bubble is ellipsoid, the supernova will first collide with outer shell at the “waist” of the bubble. 


1989 Ohno & Shibata 


1993 Shabala et al. 2010), shows 


Fig 

tact 


■ 0 , 

t ais 


a similar evolution as Simulation B, though with a more 
asymmetrical bubble. Initially, the difference between the 
5[iGand 10 pG fields is barely noticeable (extreme left of 
Fig. 0. However, as time progresses the greater strength 
of tl^ 10 pG field becomes clear (right side of left panel of 
leading to a bubble with an asymmetry at the con- 
iscontinuity of approximately 2.5:1 at the end of the 
main sequence (left side of centre panel in Fig.[^. Along the 
longer axis, the radius of the bubble is approximately 31 pc, 
which makes it slightly larger than the radius of the bubble 
in an ISM lacking a magnetic field. Along the shorter axis 
the bubble has a radius of only 12 pc. As for Simulation B, 
the outer edge of the shell (visible only as a slight distortion 
of the magnetic field lines) remains nearly spherical because 
of the pressure exerted by the compressed magnetic field. 


The strong asymmetry of the bubble influences the evo¬ 
lution of the circumstellar nebula formed during the WR 
phase (right panel of Fig.[^. In the direction perpendicular 


to the magnetic held the nebula collides with the swept-up 
ISM shell much earlier than in the direction parallel to the 
held. As a result, material is driven back toward the central 
axis where it piles ujQ 

Once the supernova explodes, the asymmetry of the 
outer bubble strongly influences its expansion. In the di¬ 
rection perpendicular to the held, the supernova remnant 
hits the outer edge of the shocked wind bubble early (left 
side of centre panel in Fig. |^. Even as it is already moving 
back toward the centre from this collision, it is still ex¬ 
panding in the direction parallel to the held (right side of 
centre panel in Fig.[^. Eventually, the returning supernova 
expansion collides with itself, along the major axis of the 
bubble. By this time most of its momentum has been lost 
and its velocity becomes subsonic (right panel of Eig. |^. 


^ This effect is exacerbated by the 2D symmetry of the simu¬ 
lations, which does not allow flow over the poles. 
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R-axis (xlO^lS cm) R-axis (xlO^lS cm) R-axis (xlO^lS cm) 

Fig. 7. Similar to Figs. and but for Simulation C with an external magnetic field of 10 [iG. This plot shows both 
the logarithm of the density in^gs] and the magnetic field lines. The scale of the central and right panels in this figure 
is 80x80 pc. As for simulation B, the swept-up shell of shocked ISM is relatively thick, which prevents the formation of 
thin-shell instabilities; and the contact discontinuity between shocked wind and swept-up ISM is strongly asymmetric. 
Although the RSG shell can form freely against the spherically symmetric wind termination shock (right side of centre 
panel), the expanding WR nebula is strongly affected by the asymmetry of the bubble (left side of right panel). 



-100 -50 0 50 100 -100 -50 0 50 100 -100 -50 0 50 100 

R-axis (xlO^lS cm) R-axis (xlO^lS cm) R-axis (xlO^lS cm) 


Fig. 8. Similar to Figs. and but for simulation C with an external magnetic field of 10 pG. Because the supernova 
expansion is spherical aim the outer boundary of the bubble is strongly ellipsoid, the supernova will first collide with the 
outer shell at the “waist” of the bubble, noticeable on the left side of the right panel. In fact, on the right side of the 
right panel, the supernova is already moving back toward the centre in the perpendicular direction, whereas it is still 
expanding in the direction parallel to the field. 


4.4. High interstellar magnetic field: 20luG 

Because of the external magnetic field strength, which, at 
20 pG, is four times as high as in simulation B, simulation D 
has from an early stage a much more asymmetric bubble 
(left panel of Fig. |^. In the direction perpendicular to the 
magnetic field, the shocked wind bubble is unable to ex¬ 
pand beyond about 7 pc. Once it is has reached this point 
the expansion stops (at least as far as the contact discon¬ 
tinuity between the hot bubble and the swept-up shell is 
concerned), because the thermal pressure in the hot bub¬ 
ble drops below the magnetic pressure of the ISM. This 
also limits the expansion of the wind termination shock, 
which reaches its maximum distance of 2.7pc (right side 
of left panel of Fig. [^. Parallel to the magnetic field, the 
expansion continues, making the shocked wind bubble in¬ 
creasingly a-spherical. As in the simulations B and C, the 


interstellar magnetic field is disturbed in a sphere around 
the shocked wind bubble, although not as outspoken as a 
swept-up shell. 


As the bubble becomes increasingly a-spherical, a pres¬ 
sure gradient starts to form along the major axis, with the 
thermal pressure in the “tips” of the bubble becoming lower 
than the ram pressure of the wind at the termination shock. 
This gives the bubbles a somewhat eye-like shape, rem¬ 
iniscent of the much smaller eye-like structure s observed 
aroun d some Asymptotic Giant Branch stars (|Gox et 


2012). To the “north” of the star, the magnetic held actu¬ 


ally squeezes off part of the bubble (centre panel of Fig.fl^, 
though it should be noted that this effect is enhanceoDy 
the 2.5-D symmetry of the model. As a result of the re¬ 
duced effective volume of the bubble, the thermal pressure 
increases temporarily, pushing the wind termination shock 
back toward the star. The resulting bubble shape is remi- 
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Fig. 9. Similar to Figs. andbut for simulation D, which has a 20 [iG interstellar magnetic field The shocked wind 
bubble is highly asymmetric from the earliest stages of the evolution and stops expanding in the direction perpendicular 
to the magnetic field after it reaches approximately 7 pc. The elongated shape of the bubble strongly influences the 
development of the WR nebula (right panel). 



Fig. 10. Similar to Figs. and but for simulation D with an external magnetic held of 20 pG. The scale of each 
panel is 120x120 pc. The strongly asymmetrical shape of the bubble causes the supernova to expand along the z-axis 
after recoiling from the outer edge. 


niscent of the “eye-l ike” shapes found for AGB stars by |van| 


Marie et al. (2014a), albeit on a much larger scale. 


When the star becomes a RSG, the ram-pressure of the 
wind decreases considerably (right side of centre panel in 
Fig.§. Since the cross-section of the bubble perpendicular 
to the magnetic held is much shorter than in the direction 
parallel to the held, the effect of change in the wind is first 
noticeable around the “waist” of the bubble. The reduced 
ram-pressure leads to a reduction in thermal pressure. This 
in turn, causes the contact discontinuity to move back to¬ 
ward the star when the thermal pressure of the shocked 
wind can no longer counterbalance the magnetic pressure 
of the ISM. The reverse happens when the star reaches the 
WR phase and the WR shell sweeps up the free-streaming 
RSG wind as well as the RSG shell at the termination 
shock. Unable to expand in the r-direction, the shells col¬ 
lide with the outer edge of the shocked wind bubble and are 
focussed along the longitudinal axis of the bubble, forming 
a jet-like structure. Afterwards, the strong WR wind causes 
an increase in thermal pressure in the shocked wind bubble. 
Once again, this is first felt at the “waist” of the bubble, 


which now expands until the magnetic held brings it to a 
halt again. The increased thermal pressure also re-opens the 
connection to the squeezed-off part (left side of left panel 
of Fig. At the end of the stellar evolution, the bubble 
is still highly a-spherical, with an asymmetry of approxi¬ 
mately 7:1. Along the axis parallel to the magnetic held, 
the radius is 50 pc, but at the “waist” the radius is only 
about 7 pc. 


Once the supernova explodes, the expanding remnant 
quickly collides with the outer shell because of the narrow 
“waist” of the bubble and recoils, flowing back to the cen¬ 
tral axis. It then expands along the length (centre panel of 
Fig. 10). Eventually, the expansion reaches the end of the 
elongated bubble and recoils toward the centre, Ailing up 
the cavity (right panel of Fig. 10). 
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Fig. 11. Similar to Figs. 1^ and § but for simulation E, which combines a 5 [JiG magnetic field with a warm, low 
density medium. The bubme expands very quickly and becomes much larger than for the previous simulations. (Note 
the different scale: the left panel is 160 by 160 pc, the centre and right panel are 290 by 290 pc.) This allows the red 
supergiant and Wolf-Rayet shells to form unimpaired by the aspherical shape of the bubble. 



Fig. 12. Similar to Figs. 1^ and i but for simulation E, which combines a 5 pG magnetic field with a warm, low 
density medium. All panels are 290 by 290 pc. The impact of the supernova remnant with the contact discontinuity 
causes the local gas to be compressed, which leads to the formation of instabilities (right panel). 


5. Bubbles in a warm ISM 


5.1. The influence of a warm, low density ISM 


Simulation E (Figs.|ll||12|) combines a warm IMS (10 000 K) 


with low density (10 ^ g cm“^) and a 5pG, magnetic 

field. Because both the ISM thermal pressure and the ISM 
inertia are low, owing to the low density, the bubble can 
expand very quickly (left panel of Fig. 11), remaining su¬ 
personic for the early part of the main sequence phase, de¬ 
spite the higher sound speed in the ISM. This leaves the 
magnetic field as the dominant force in the ISM, though 
it is insufficient to stop the expansion. The end result is 
a very large, ellipsoid bubble. By the end of the main se¬ 
quence (left side, centre panel of Fig. 11), the bubble has 
reached a size of approximately 270 by 38 parsec. Because 
of the extremely large size of the bubble, the development 
of the internal shells (RSG shell and Wolf-Rayet ring neb¬ 
ula) proceeds unimpaired, with the Wolf-Rayet ring nebula 
retaining its spherical shape and breaking up before it col¬ 
lides with the outer edge of the main sequence bubble (right 
panel of Fig. [n]). 

The supernova expansion also proceeds largely uninhib¬ 
ited (left and centre panels of Fig. 12). However, unlike the 
previous simulations (A-D), the expansion of the supernova 
remnant actually leaves a visible imprint on the interface 


with the ISM, creating a thin, localized shell and some de¬ 
formation of the bubble (right panel of Fig. 12). This is the 
result of the low ISM density, which allows the contact dis¬ 
continuity to accelerate outward under the force excerted 
by the collision with the supernova remnant, compressing 
the interestellar gas into a thin shell, despite the magnetic 
field pressure. This also leads to the formation of Rayleigh- 
Taylor instabilities at the contact discontinuity (right side 
of right panel of Fig. 12). 


5.2. The influence of a warm, dense ISM 

For simulation F, we repeat simulation B but with an ISM 
temperature of 10 000 K. This type of combination (both 
warm and high density) is expected to reside in Hu regions 
in the close vicinity of massive stellar clusters. It serves 
to demonstrate how the interstellar magnetic field and ISM 
thermal pressure combined can influence the shape and size 
of the circumstellar bubble, leading to a radically different 
bubble evolution. 

Under these conditions the thermal pressure in the ISM 
(2.6 X 10“^^ dyne/cm^) is more than an order of magnitude 
higher than the magnetic field pressure. The effect of the 
two pressures combined leads to a highly a-typcial bubble 
development as shown in Figs. pMjfM) Due to the high ther- 
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Fig. 13. Development of a circumstellar bubble during the main sequence in a high density, warm ISM (simulation F) 
with a 5pG interstellar magnetic field. From left to right this figure shows density and magnetic field lines at 1,2, 4.3, 
and 4.5 Myr. The bubble is much smaller compared to the other SpG simulations (B & E). The combined thermal 
and magnetic pressure in the ISM is strong enough to reduce the size of the bubble considerably. By the end of the 
main sequence (left side of right panel), the bubble has obtained a somewhat “eye-like” shape. At the onset of the red 
supergiant phase (right side of right panel), the hot, shocked wind bubble collapses under the combined forces of the 
interstellar magnetic field and the thermal pressure in the ISM. 



R-axis (xlO'^lS cm) 
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R-axis (xlO^lS cm) 


Fig. 14. The shape of the bubble for simulation F after 4.45, 4.5, 4.55 and 4.6 Myr. During the red supergiant phase 
(left panel), the hot, shocked wind bubble is destroyed completely by the combined forces of the interstellar magnetic 
field and the thermal pressure in the ISM. The result is a purely radiative collision between the red supergiant wind and 
the outer shell (left side of left panel.) Once the start reaches the Wolf-Rayet stage (rightside, left panel), the powerful 
wind quickly sweeps up the red supergiant wind, forming a new hot bubble that starts to expand into the ISM (right 
panel). 


mal pressure in the ISM the expansion is slow, being both 
sub-sonic and sub-Alfvenic, which reduces the internal vol¬ 
ume of the bubble. This in turn increases the internal den¬ 
sity, making radiative cooling more efficient, which reduces 
the internal pressure of the bubble, leading to a similar 
process as in simulation D: the outer tips of the expanding 


bubble are squeezed off by a combination of thermal and 
magnetic pressure from the ISM, resulting in a somewhat 
“eye-shaped” bubble ( left side of right panel o f Fig, 


et al. 


2012 
that t 


similar to the ones found aroun d some AGB stars ( Gox 

|van Marie et al.||2014a), though it should be noted 
;his effect may be an artifact of the 2D nature of our 
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Fig. 15. Similar to Figs. ^^^ 


10 


and 


12 


but for simulation F. This figure shows the expansion of a supernova inside a 


wind bubble that is constrained by both high thermal pressure and magnetic field pressure. The supernova is effectively 
trapped inside a relatively small bubble and, as in all previous models, fails to break out of the shell. 



R-axis (xlO^ia cm) 



R-axis (xl0''18 cm) 



R-axis (xlO''ia cm) 


simulations (see Section 6.6). Afterwards the bubble starts 
to expand once more because the energy from the wind is 
now contained in a relatively small volume. However, be¬ 
fore it can regain its former size, the star reaches the end of 
the main sequence and becomes a red supergiant (right side 
of right panel of Fig. 13). The resulting decrease in wind 


ram pressure stops the expansion and causes the bubble to 
collapse completely under the outside forces ( right side or 
right panel of Fig, and left panel of Fig.[^. The result is 
a fundamental change in the interaction between the wind 
and the ISM. Where before the wind termination shock was 
adiabatic, it now becomes completely radiative with the 
free-streaming red supergiant wind impacting directly onto 
the ISM. By the time the star evolves from red supergiant 
to Wolf-Rayet star the hot, shocked wind bubble has com¬ 
pletely disappeared ( left panel of Fig. pM|). The Wolf-rayet 
wind forms a ring nebula around the star (right side of left 
panel of Fig. 14), which, driven by the hot, shocked Wolf- 


Rayet wind, quickly expands into the red supergiant wind, 
forming a new, hot shocked wind bubble, which continues 
to expand outward, effectively recreating the original main 
sequence bubble. The expansion is, once again, somewhat 
a-spherical, owing to the magnetic field pressure as well as 
the fact that the new expanding shell encounters piled-up 
mass along the z-axis. By this stage the bubble closely re¬ 
sembles the results that van Marie et al. (2014a) found for 
a planetary nebula expansion inside an AGB wind that had 
been constrained by an interstellar magnetic field. 

By the time the star explodes as a supernova (Fig. 
the bubble is still small compared to the other bubbles .^s 
a result, the supernova remnant quickly reaches the outer 
shell and bounces back. Because the bubble is not extremely 
aspherical, the interaction between the supernova remnant 
and the outer shell occurs almost simulteaneously along 
the entire shell and avoids the strongly asymmetrical shape 
that occured in simulation D. 


6. Influence of interstellar magnetic fields on the 
properties of circumstellar bubbles 

6.1. General shape and size of the circumstellar bubble 


As predicted analytically ( Heiligman|198Q ) and shown pre- 
viously in numerical models by Tomisaka| (|199Q[ |1992|) and 
Ferriere et al. (1991), the magnetic held limits the expan¬ 


sion of the bubble in the direction perpendicular to the 
field, partially because of the magnetic tension force, which 
counteracts the bending of the field lines by the expand¬ 
ing bubble, and partially by the magnetic pressure, which 
counteracts the field lines being compressed. 

A direct comparison between the four bubbles can be 
made through cross-sections of the bubbles as shown in 
fig. which show the density (left panel) and magnetic 
field^rength (right panel) along 1-D cuts perpendicular to 
the field (continuous lines) and at a 45° angle (dashed lines) 
for all four simulations at the end of stellar evolution. In all 
cases, the shell is characterized by a thin, highly compressed 
feature at the contact discontinuity, where the shocked 
wind directly sweeps up the surrounding medium and an 
extended, only slightly compressed outer shell, where the 
magnetic field, compressed by the expanding bubble pushes 
the gas ahead. For increasing magnetic field strength, the 
amount of compression perpendicular to the field in the 
outer shell decreases, until, at a field strength of 20 pG, one 
can barely speak of an outer shell. This is the consequence 
of the increase in the Alfven speed, which is defined as, 

B 

~ 74 ^' 

Even for a 5 pG field in an ISM with a density of 
10“22-5 gcm“^, this gives us an Alfven speed of 2.5kms“^, 
which is significantly higher than the adiabatic sound speed 
in a cold ISM (^ lkms“^ at 100K). Hence, the strength 
of the shock is determined by the Alfven speed, rather 
than the sound speed. A higher Alfven speed reduces the 
strength of the shock and therefore the compression rate. 
At 20 pG, the Alfven speed is approximately 10kms“^. At 
the same time the total expansion of the bubble, perpendic¬ 
ular to the field, at 20 pG is only 7 pc in 4.77 million years: 
1.4kms“^. Therefore, the expansion perpendicular to the 
field has become a sub-Alfvenic wave, rather than a shock. 
The same effect, though to a lesser extent, is noticeable at 
a 45° angle to the field. 

The behaviour of the bubbles in simulations B-D gen¬ 
erally follows the predictions based on Fig. For those 
circumstances where the ram pressure of the stellar wind is 
always higher than the magnetic pressure of the interstellar 
magnetic field, as is the case for the 5 pG model, the bubble 
becomes ovoid, rather than spherical owing to the magnetic 
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Fig. 16. These figures show the gas density (left) and mametic field strength along cuts through the 2D simulation results 
at the end of stellar evolution (far-right panels of Figs. [8| and [Tq|). Each cut starts at the central star and moves 
either along the R-axis (dashed lines) or at a 45° angle (continuous lines) between the R- and Z-axes. The compression 
of the outer shell decreases with increasing field strength, Alfven-waves weakening the shock. 


tension force that seeks to reduce the curvature of the field 
lines, but the expansion continues in all directions. If the 
magnetic field pressure exceeds the ram pressure, expansion 
perpendicular to the field lines ceases altogether, resulting 
in a tube-like bubble such as we find for the 20 pG model. 
The 10 pG model falls in between these two, but generally 
follows the 5 pG model, because the ram pressure of the 
wind only briefly drops below the magnetic field pressure. 

As soon as Pram < Pb expansion ceases in the direction 
perpendicular to the magnetic field, at least for the inner 
boundary of th e swept-up ISM shell. T his d id not occur in 


the m odels by Tomisaka (1990 1992) and Ferriere et al. 
(1991), who did not investigate magnetic fields of sufficient 
strength, and in our simulations happens only for the model 
with PisM = 20 pG. The outer limit of the shell continues 
to expand, driven by the magnetic pressure, which reduces 
the local density inside the shell. Interestingly, owing to 
the extremely wide “shell”, the field lines are affected to a 
distance of about 50 pc, despite the fact that the tube-like 
shocked wind bubble is only about 7 pc across. As calcu¬ 
lated in Sect. for our main sequence wind parameters 
the radius for the termination shock should be 2.6 pc if the 
interstellar magnetic field strength is 20 pG. This is well in 
agreement with our numerical results, which put the wind 
termination shock at a maximum of 2.7pc (right side of 
the left panel of Fig. [^. In theory, the Wolf-Rayet wind is 
strong enough for the expansion perpendicular to the field 
lines to resume, even for the 20 pG magnetic field. We see 
some sign of this, with the “waist” of the bubble expanding. 
However, since the Wolf-Rayet phase is relatively short, it 
lacks the time to change the general tube-like appearance 
of the bubble. 


As was shown for super bubbles by Tomisaka 
Tdi 


1992^ and Ferriere et al. 


1990 

19911), the total volume of the in¬ 


dividual bubbles is greatly reduced by the presence of the 


magnetic ISM field. The absolute size of the major axis 
of the bubble does not change much with the magnetic 
field strength for those simulations where the magnetic field 
only slows the expansion (simulations B and G). However, 
should the field stop completely the expansion in the per¬ 
pendicular direction, as is the case for simulation D, then 
the excess energy is used to expand further in the direc¬ 
tion parallel to the field, increasing the asymmetry of the 
bubble. 

In the non-magnetic case, the final bubble is spheri¬ 
cal and has a radius of approximately 30 pc. This gives 
the shocked wind bubble a final volume of more than 
110 000 pc^. For the 20 pG interstellar magnetic field, the fi¬ 
nal bubble is almost cylindrical, with dimensions r = 7 pc 
and z = 100pc, for a total volume of about 15 000pc^. 
As a result, the density of the shocked wind is eight times 
higher for the magnetic model. This makes the radiative 
cooling in the shocked wind 8^ = 64 times as effective, 
further reducing the thermal pressure in the shocked wind. 
(The length-width ratio for the bubbles inside an interstel¬ 
lar magnetic field may have been exaggerated by the fact 
that our 2D simulations enforce axi-symmetry, which in 3D 
could be broken due to local effects.) 

In the low density, warm, ISM (E), the magnetic field 
dominates the force on the outside of the expanding bub¬ 
ble (owing to a lack of either thermal pressure or inertia). 
As a result the bubble becomes far more asypherical than 
when a similar magnetic field is combined with high density 
gas (B). Furthermore, the low density allows the bubble to 
grow to a much larger size. In other respects, the bubble 
evolution still follows the same pattern as in simulations B 
and G. 

For the special case where a warm, high density ISM is 
combined with an interstellar magnetic field (simulation F), 
the evolution of the circumstellar bubble deviates radically 
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from the established norm; particularly during the red su¬ 
pergiant phase, when the hot, shocked wind bubble disap¬ 
pears completely, to be replaced with a purely radiative 
shock between the red supergiant wind and the ISM. 


6.2. From circumstellar bubbles to super bubbles 

Although we limit ourselves to the circumstellar medium 
of one star, in reality massive stars are formed in clusters. 
This changes the evolution of the circumstellar medium, 
because the wind-blown bubbles of the individual stars en¬ 
counter one another as they expa nd and merge to f orm a 


Besides bubble dynamics and morphological aspects, 
the interplay of bubbles and their turbulent environment 
i s essential to understand the propagation of cosmic rays 


and Krause et al. (2013a|b 

). As demonstrated both in this 

paper, as well as in van IV 

arle et al. 

(2012), each circum- 


stellar bubble can effectively contain the expansion of the 
supernova, thus even when the stars explode no interaction 
occurs, unless the bubbles hav e already merged durin g the 
wind-driven expansi on phase. Tomisaka ( 199Q[ 1992) and 
Ferriere et al. (1991) explored the influence of tne interstel- 


lar held on the outer shape of a super bubble, but could not 
model the mergers, because they approximated the cluster 
as a point source of kinetic energy. 

The presence of a strong interstellar magnetic field could 
potentially prevent such mergers, depending on the num¬ 
ber of stars, their distribution over the local space and the 
strength of the magnetic field. Along the field lines, the in¬ 
dividual bubbles will merge easily. However, perpendicular 
to the field lines, the bubbles will have to overcome the in¬ 
terstellar field, the strength of which is enhanced because 
the field is compressed from multiple directions at the same 
time. The assumption of a large scale, uniform magnetic 
field would no longer be correct under these circumstances. 
However, such a field is not necessary to influence the shape 
of the circumstellar bubbles. In a dense cluster, the local 
ISM density at the start of star formation will most likely 
be much higher than 10cm“^. At the same time, the mag¬ 
netic fiel d is likely to be at least an order of magnitude 
stronger (Crutcher 1999). As long as the field is uniform 
on the same scale as the individual bubbles (a few parsec 
for ISM densities of molecular clouds) the expansion of the 
bubbles can be contained. The final outcome would depend 
on the number and density of the star cluster. If the indi¬ 
vidual bubbles are kept from merging, this would result in 
a filamentary structure, not only in terms of density, but 
also in terms of metallicity, because the enriched supernova 
material would effectively be trapped inside the elongated 
wind bubbles and only merge with the ISM when the bub¬ 
bles collapsed due to a lack of internal pressure after the 
deaths of their progenitor stars. This complex investigation 
will be carried out in further works. 

The general shape of superbubbles can serve as a test 
for the influence of the interstellar magnetic held on the 
expansion of wind-driven bubble s, even if the held does 
not prevent mergers. According to Heiles (1979) super bub- 
bles tend to be elongated along the galactic plane. [Ferriere 


[Bykov 


2001 


et al. (1991) argued that this indicated that magnetic fields, 
rather than density stratification dominated the shape of 
super bubbles, because the magnetic held orientation in the 
galaxy tends to be aligned with the disk, whereas density 
stratification would increase the expansion perpendicular to 
the disk. A similar alignment was observed by jWeidmannj 


Ferrand & Marcowith 2010). Cosmic rays, 
epending on the way they are transported inside these 
structures, may also impact the bubble dynamics; a pos¬ 
sibility that has been largely overlooked until now. These 
complex aspects deserve a full 3D approach including the 
impact of energetic particles. Software capable of handling 
such a model is currently under development. 


6.3. Effects on the shape of shells inside the bubble 

Although the held cannot affect the events inside the bub¬ 
ble, such as the formation of circumstellar nebulae and su¬ 
pernova remnants, directly, it can affect them indirectly, 
because these events are constricted by the size and shape 
of the bubble. 

For the weak (SpG) and intermediate (10 pG) fields the 
shells formed inside the bubble as a result of changing wind 
parameters are largely unaffected by the presence of the 
interstellar magnetic held. The cross-section of the bubble, 
even along the minor axis, is large enough to allow the shells 
to form and evolve as though the whole bubble were spher¬ 
ical. For the strong (20 pG) magnetic held the situation is 
different. In this case the bubble has become so constricted 
by the magnetic held, that the evolution of the WR shell is 
impacted at an early stage and forced into an asymmetric 
shape that resembles a jet, even though it was originally 
formed as a spherical shell. The same holds true for the 
shells formed by the expansion of supernova remnants. 

The influence of the interstellar magnetic held on the 
shape of supernova remnants can potentially serve as a 
means of testing the model observationally. The interstellar 
magnetic held cannot influen ce the shape of an expanding 
supernova blastwave directly (Manchester 1987). Therefore, 
any correlation between the shape of supernova remnants 
and the morphology of the interstellar magnetic held has 
to be caused indirectly: The magnetic held shapes the wind 
cavity into which the supernova explodes and thereby in¬ 
fluences the shape of the supernova remnant indirectly. 
That the wind-cavity influences the shape of the supernova 
re mnant was first demonstrated throu gh 2D simulations 
by Rozyczka et al. (1993). Arnal ( 1992|) proposed ambient 
magnetic fields as an explanation for barrel shaped cavi¬ 
ties around Wolf-Rayet stars, which in turn cou ld be the 
source of asymmetry in supernova remnants and Gaensler 
(1998) showed that galactic supernova remnants tend to 
be aligned with the galactic disk, which seems to confirm 
that the interstellar magnetic held affects their expansion. 
A drawback of this test is that an asymmetric cavity around 
a Wolf-Rayet star (and as a result an asymmetric supernova 
remnant) can in theory be formed through the motion of 
its progenitor star relative to the surrounding ISM. If the 
motion is supersonic, it causes the stellar wind to form a 
bowshock structure, rather than a normal bubble. Although 
no proven case of a bow shock around a Wolf-Rayet star 
has been observed, it can theoretically accou nt for the ex¬ 
istence of asymmetric supernova remnants (Gvaramadze 
2006; Meyer et al.||2Q15 ). 


& Dfaz (2008) for planetary nebulae. 


In simulation E the a-spherical shape of the bubble is 
largely irrelevant to the development of the internal shells, 
on account of its shear size. Even the Wolf-Rayet ring neb- 
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ula has broken up and largely dissipated before it collides 
with the outer edge of the bubble. 

Again, simulation F, with its unique bubble evolution 
presents a special case. Because the main sequence wind 
bubble disappears completely, there is no separate red su¬ 
pergiant shell. Instead, the red supergiant wind impacts 
directly on the ISM. A Wolf-Rayet ring nebula is formed 
inside the red supergiant wind, but does not survive the 
collision with the outer shell. 


6.4. Lack of instabilities 

In all three simulations including a magnetic field (B, C 
and D) the magnetic field effectively prevents the forma¬ 
tion of instabilities in the outer shell. The magnetic pres¬ 
sure counteracts the compression of the shell, preventing 
the formation of thin-shell instabilities . In theory, Rayleigh- 


Taylor instabilities could still form ( Breitschwerdt et al 
20001 iStone & Gardiner| |2007|) and we see some sign of 


the onset of such instabilities along the contact disconti¬ 
nuity between shocked wind and shocked ISM (see centre 
panels of Figs. [^and[^. Rayleigh-Taylor instabilities are 
also evident in Simulation E when the interstellar gas at 
the contact discontinuity is temporarily compressed by the 
impact of the supernova remnant (right panel of Fig. 12). 
However, further growth of such instabilities would involve 
motion of the gas perpendicular to the field lines and the 
forces creating the instability are insufficient to overcome 
the magnetic field. This side-effect may be of interest in 
other circumstances as well and can, for example, explain 
the absence of instabilities in the Herschel observations of 


the bow-shock of a-Orionis (Decin et al. 2012 


van Marie 


et al. 2014b|). When studying the effect of the interstellar 
magnetic field on the instabilities, one should keep in mind, 
that our models are limited by their 2.5-D nature. The in¬ 
teraction between a magnetic field and hydro dynamical in¬ 
stabi l ities is an inherently 3D problem (Breitschwerdt et al. 
2000 Stone & GardinerJ|2007) and should be further inves- 


ti gatedTetore conclusi o ns can be drawn (see the discussions 
in Arthur et al. (2011); Kim & Kim ( 2014| )). 


6.5. Implications for gamma-ray burst progenitors 

Because the magnetic field keeps the bubble confined, the 
thermal pressure of the shocked wind remains higher, which 
places location of the wind termination shock, the balance 
between thermal pressure and ram-pressure, closer to the 
star, though not as high as predicted by the analytical ap¬ 
proximation (See Sect. |^. This can be of interest in the 
case of possible long gamma-ray burst (GRB) progenitors. 
The afterglows of many long GRBs show evidence of a con¬ 
stant density, rather than wind-like environment, despite 
the fact that these bursts are thought t o originate from 


massive stars (e.g. Ghevalier et al. 2004). This has been 


explained as being the shocked wind, rather than the free- 
streaming wind. However, this would indicate that the wind 
termination shock is very close (^0.1 pc) to the star. This 
has bee n explained in several ways, including a h igh ISM 
density ( Eldridge et al.|20Q6| [van Marie et al.|20Q6|), and/or 
a high ISM temperature dGhevalier et al.||2004[ |van Marie 


et al. 2006). However, it is possible that interstellar mag- 
netic helds play a role here as well. Although the mag¬ 
netic fields applied in this paper do not reduce the radius 


of the wind termination shock as much as is required to ex¬ 
plain the GRB afterglow observations, higher mag netic field 
strengths h ave been observed in molecular clouds (Crutcher 


1999| 2012|). These fields, combined with other factors such 
as high ISM density and temperature, and a reduced mass 
loss rate at low metallicity, could, in theory, reduce the ter¬ 
mination shock radius to less than 0.1 pc. 


6.6. Limitations of the model 

We are aware that the ISM parameter range explored in 
this work is limited. The environment of massive stars is 
very complex, inhomogeneous, and intermittent. The struc¬ 
ture of the parent molecular cloud includes media of dif¬ 
ferent properties from diffuse to dense cores. The molecu¬ 
lar cloiM^_axe_^nown__to_^ turbu lent (see a recent review 
by Hennebelle & Ealgarone (2012)), pervaded by magnetic 
fields, the stre ngth of which ca n exceed 20 pG by an order 
of magnitude (Crutcher 2012). The turbulent motions in 
the interstellar medium are known to have an imp act on 
the shape and the dynamics of t he stellar bubble (Korpi 


et al. 


1999 Silich fc Eranco|p~999 ). Hn regions also harbor 


turbulentmotions sh owin^articular statistical properties 
(Medina et al. 2014); turbulent motions associated with 


various processes as: internal globules evaporation and/or 
interaction with stellar winds. 

More realistic models of the ISM into which the circum¬ 
stellar bubbles of massive stars are blown have to include 
the combination of magnetic field effects in turbulent molec¬ 
ular clouds and/or in Hu regions. Such models deserve a full 
3D treatment. Eurthermore, before such comprehensive 3D 
models are attempted, it is necessary to first explore how 
each of the many interacting forces in the ISM affects the 
evolution of a circumstellar bubble. To this end we have 
deliberately kept the gas properties of the ISM a constant 
in order to isolate the influence of the large scale magnetic 
field on the evolution of the circumstellar bubbles. 

The 2D nature of our models, while greatly reducing 
computation time, forms a second limitation. Eor purely 
hydrodynamic interactions van Marie & Keppens (2012) 
showed that instabilities in Wolf-Rayet ring nebulae tend 
to grow more quickly in 3D models than in 2D, though 
the difference tends to decrease over time. As to the ef¬ 
fect of 2D symmetry on the magneto-hydro dynamical in¬ 
teraction between the bubble and the interstellar magnetic 


field, this is a more complicated question, [van Marie et al. 


(2014a|) showed that a 3D model of the interaction between 
an AGB wind and the interstellar magnetic field yields the 
same qualitative and quantitative result as a 2D model. 
However, that model lacked some of the features that can 
occur in the models presented here: in particular the phe¬ 
nomenon observed in simulations D and E, where part of 
the bubble is separated from the rest because of the mag¬ 
netic field pressure. This effect is particularly visible in sim¬ 
ulation E during the late main sequence (Eig. [T7|. As the 
bubble expands the thermal pressure in the regions far from 
the star decreases to the point where it becomes less than 
the combined magnetic and thermal pressure of the ISM. 
As a result these regions are cut off from the rest of the 
bubble and disappear since they no longer have a stellar 
wind to supply them with internal energy. In a 3D model, 
which has an additional degree of freedom for expansion, 
this effect may be reduced, or even disappear completely. 
If it occured at all, it would not show the axial symmetry 
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Fig. 17. The shape of the bubble for simulation F after 3, 3.5, 3.75, and 4Myr. The tips of the bubble are squeezed off 
by the magnetic field. This effect is enhanced by the 2D nature of our simulations and may not actually occur in 3D. 


that occurs in 2D. The “eye-like” shape of the bubble after 
the tips have been reduced is not a 2D artefa ct and has 
been observe d in 3D simulations of AGB winds Fvun Mar lei 
et al.l (|2014a|) 


So far we have ignored the existence of stellar magnetic 
fields. Simulations of the circumstellar medium of low mas 


stars (e.g. Rozyczka & Franco 1996 Garcfa-Segura et al. 
1999[ ) show that toroidal stellar magnetic helds of AGB 


and post-AGB stars can create an extraordinary variation 
in the shapes of planetary nebulae by changing the direc¬ 
tion of the stellar wind, which for these kind of stars tend 
to be weak. Simulations of the interaction between the stel¬ 
lar wind and magnetic field of high mass stars tend to fo¬ 
cus on dipole fields in a region close to the star (primar¬ 
ily the wind acceleration zone) rath er than the extended 
wind cavity. Such simulations (e.g. Ud-Doul a & Owocki 
' ^ lUd-Doula et al. 


2002 Townsend et al. 20071 


show that while a dipole stellar magnetic held can 


'i 


| |2UU8[ 
ild can 


200 ^ 

influ' 


ence the morphology of the stellar wind close to the star, 
at longer distances the stellar wind, which is faster and far 
more powerful than the winds of low mass stars, starts to 
dominate, forcing the field lines to align with the expand¬ 
ing wind to form a magnetic monopole. This allows the 
wind to expand spherically as though the magnetic field 
no longer exists, leaving the magnetic field largely irrele¬ 
vant to the morphology of the circumstellar shells. Exactly 
how such a magnetic field would interact with the inter¬ 
stellar magnetic field is difficult to predict, but it should 
be noted that the distance from the star would greatly 
reduce the s tellar magnetic field strength. For example, 
Wade et al. (2012) estimate the surface dipole strength of 
NGC 1624-2, the strongest magnetic field ever measured 
for an 0-star, to be at least 20 kG. Assuming that the field 
strength decreases with the distance squared (as would be 
the case for a mono pole) and a stellar radius of about 11 Rq 
( Wade et al.| 2012] ) this would reduce the field to approx- 
imately 1.2 x 10“^ pG at a distance of one parsec, sig¬ 
nificantly weaker than even the weak interstellar fields in 
the galactic halo. For more representative 0-star magnetic 


field strengths of 10-1000 G Walder et al. (2012 and refer¬ 
ences therein) the effective field strength at the interaction 
point with the ISM would be even weaker. Note that if the 
stellar magnetic field were to retain its dipole (or higher or¬ 
der) configuration despite the stellar wind, its field strength 
would decrease even faster. Gonsidering that wind-cavities 
tend to be on a multi-parsec scale, it seems reasonable to 
assume that the interaction between the stellar and inter¬ 
stellar magnetic fields would only be relevant for very young 
stellar objects and become irrelevant once the wind cavity is 
formed. However, these simulations only involved hot stars 
such as O- and B-stars. During the red supergiant phase 
the ram pressure of the wind is much smaller and the wind 
velocity is lower, reminiscent of an AGB star. For such as 
star a toroidal magnetic field, if present, may influence the 
direction of the wind in a manner similar to that predicted 
for AGB stars. However, the scale of the circumstellar cav¬ 
ity (at least an order of magnitude large than those of low 
mass stars) would still limit any direct interaction between 
the stellar magnetic held and the ISM held. 


7. Conclusion 

The shape, size, and evolution of circumstellar bubbles are 
influenced by many factors, such as the properties of the 
stellar wind, the motion of the star relative to its surround¬ 
ing, the density and temperature of the interstellar gas, and 
the strength and shape of the interstellar magnetic held. In 
this paper we have focussed our investigation on the influ¬ 
ence of the strength of the interstellar magnetic held on the 
shape and size of the circumstellar bubble around a massive 
star. 

We have limited our calculations to regular magnetic 
fields of different strengths that are characteristic of dif¬ 
ferent regions in the galaxy: from low magnetic found in 
the disk to higher values found in the galactic bulge. The 
simulations are 3D axi-symmetric and permit a fast explo¬ 
ration of the parameter space. Although simple with re¬ 
spect to any realistic interstellar medium, the simulation 
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setups allowed us to derive some important conclusions. 
As expected, the interstellar magnetic field constrains the 
expansion of a circumstellar bubble in the direction perpen¬ 
dicular to the field. If the magnetic field is strong enough 
to generate a pressure that exceeds the ram-pressure of the 
stellar wind at the wind termination shock, the expansion 
in the perpendicular direction stops completely. This leads 
to the formation of a strongly asymmetric bubble, the size 
of which is so constrained that it affects the evolution of 
temporary circumstellar features such as Wolf-Rayet ring 
nebulae. Even for the powerful wind of an 0-type massive 
star, galactic magnetic fields can be strong enough to do 


so. 

In the future we will extend these simulations to 3D in 
order to include random fluctuations in the magnetic held. 
These fluctuation s can be of similar strength to the uniform 
held component ( Beck|[2QQ9 ). Such 3D models will also al¬ 
low us to further investigate the influence of the magnetic 
held on the local instabilities in the swept-up shell. These 
3D simulations will be complemented by a setup that in¬ 
cludes the detailed structure of the star bubble environ¬ 
ment that comprises the parental molecular clouds and Hn 
regions. 
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Fig. A.l. Animation A.wind, shows the evolution of the 
circumstellar bubble according to simulation A (no inter¬ 
stellar magnetic field) for the pre-supernova phase. 

Fig. A.2. Animation A_sn shows the expansion of the su¬ 
pernova remnant for simulation A. 

Fig. A.3. Animation B_wind shows the pre-supernova evo¬ 
lution of the circumstellar bubble according to simulation B 
(^iSM = 5pG). 

Fig. A.4. Animation B_sn shows the expansion of the su¬ 
pernova remnant for simulation B. 

Appendix A: Animations of the evolution of 
circumstellar bubbles 

This appendix contains the results of all our simulations in 
animation form. Keep in mind that the frame rate differs 
between the stellar evolution and supernova phases because 
of the different timescales involved. 


Fig. A.7. Animation D.wind shows the pre-supernova evo¬ 
lution of the circumstellar bubble according to simulation D 
(^iSM = 20 pG). 

Fig. A. 8 . Animation D_sn shows the expansion of the su¬ 
pernova remnant for simulation D. 

Fig. A.9. Animation E_wind shows the pre-supernova evo¬ 
lution of the circumstellar bubble according to simulation E 
(^ISM = 5pG, low density warm ISM). 

Fig. A. 10. Animation E_sn shows the expansion of the su¬ 
pernova remnant for simulation E. 

Fig. A.11. Animation F_wind shows the pre-supernova 
evolution of the circumstellar bubble according to simu¬ 
lation E {Bism = 5pG, high density warm ISM). 

Fig. A.12. Animation F_sn shows the expansion of the su¬ 
pernova remnant for simulation E. 


Fig. A.5. Animation G_wind shows the pre-supernova evo¬ 
lution of the circumstellar bubble according to simulation G 
(^iSM = 10 pG). 


Fig. A.6. Animation G_sn shows the expansion of the su¬ 
pernova remnant for simulation G. 



